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1. Introduction 



Diffractive heavy vector meson production was initially evaluated to leading-log accu- 
racy in energy [1], while in [2] all vector mesons were treated to the leading-log accuracy 
in In Q 2 and In 1/x. It is now well understood in QCD, to leading twist accuracy thanks 
to a generalization of the QCD factorization theorem to hard exclusive processes [3]. 
So far, for heavy quarkonium production, this knowledge has only been used to con- 
front photo- and electroproduction of J/ip [4-6]. Unfortunately, the light mass of the 
charm quark leads to many theoretical uncertainties in the final result (relativistic ef- 
fects in the wavefunction, scale uncertainties, unitarity corrections...). By studying 
T-production, which has recently been observed in photoproduction for the first time 
at HERA [7], within the same framework, one may hope to pin down some of these 
uncertainties. 

The current paper is an extension of the work presented by Frankfurt et al [5]. One 
of the main findings was that the use of realistic light-cone wavefunctions for the heavy 
vector mesons, ipv{ z , k t ), with significant average kt, leads to an overall suppression of 
the cross-section relative to the static, ipv( z ik t ) = S(z — 1/2) S^(k t ), and Gaussian, 
i[)y(z,kt) = A5(z — 1/2) exp(— akf/m 2 .), forms used in [6]. The hybrid wavefunctions 
used were designed to interpolate hard QCD behaviour at small transverse distances 
(normalised to the decay width into leptons) with quarkonium models at large trans- 
verse distances. In the limit m 2 — > oo (but m 2 <C W 2 ) such Fermi- motion effects of 
the quarks can be substantiated in QCD, because non-quark degrees of freedom will 
be suppressed by the powers of m 2 . Whether such an approximation is applicable to 
the production of states in the ^-family is an open question. The fct-suppression was 
tempered by an enhancement due to a rescaling (beyond leading-log(<5 2 )) of the gluon 
density to higher scales, to reflect more accurately the typical transverse size of the 
scattering dipoles. 

Here, we present predictions for diffractive (also called exclusive) photoproduction 
of the T-family using similar hybrid wavefunctions. One of our main results is that 
the non-diagonal (hereafter skewed) kinematics lead to a significant enhancement of 
the cross section. We present an estimate of the size of this effect, to the leading-log 
accuracy with which the skewed splitting functions are known. At HERA energies 
T(ns)-states are produced at large effective scales (around 40, 60 and 75 GeV 2 for 
n = 1,2,3, respectively) and at relatively high-x (between 0.001 and 0.02). It follows 
that the real part of the amplitude is large and we calculate it using a fit to the energy 
dependence of the imaginary part and dispersion relations. Taking both effects into 
account leads to cross-sections which concur with the measured ones and rise very 
steeply with W 2 , the 7P centre-of-mass energy (approximately W^ a85 )). 

This note is organised as follows. Firstly, we very briefly recap the relevant equations 
in [5]. Secondly, we discuss the rescaling procedure in detail. We then explain how 
the skewedness is calculated and implemented. Next, we give an explanation of the 
calculation of the real part of the amplitude. Finally, we present and discuss the 
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Figure 1: Photoproduction of Heavy Vector mesons 
resultant cross-sections, and conclude. 



2. Formulae for cross sections 



Hard exclusive processes factorize in QCD [3] and at small-x are driven by the exchange 
of two gluons in the t-channel. One such process is the diffractive electroproduction 
of heavy vector mesons, in which the exchanged gluons are connected to the gg-pair 
fluctuation of the photon (virtuality q 2 = —Q 2 ) in the four possible ways (one of which 
is shown in fig.(l)) and convoluted with ipv(z, h)- The forward differential cross-section 
for a vector meson, mass My, containing a (current) quark of mass m, can be written 
as the product of an asymptotic expression and a finite Q 2 correction, C(Q 2 ) (see [5] 
for more details) : 
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where Qljj-(Q 2 ),T V , e are the effective transverse scale, the leptonic decay width and 
photon's polarisation respectively; (3 is the ratio of real to imaginary parts of the 
amplitude. The /^-suppression factor T, R and r\y contain integrals involving the light- 
cone wavefunctions of the photon and vector meson. Hybrid wavefunctions are used 
which are designed to have the characteristic z(l — z) behaviour at very small transverse 
distances expected from QCD (this implies rj v — 3). This justifies the use of the QCD 
running mass and hence the correction for this in eq.(2.2). Overall the use of the hybrid 
wavefunctions, in place of the pure non-relativistic wavefunctions, does not lead to 
significant effects in the cross sections. In this paper we use the particular model of the 
quarkonium wavefunction due to Buchmueller and Tye [8]. They used a QCD-inspired 
potential, with a Coulomb piece corresponding to t-channel gluon exchange (oc 1/r) 
and a confinement piece (oc r). The analysis produced a value of mi, = 4.88 GeV for 
the quark mass. The significant QCD radiative corrections to the amplitude of heavy 
quarkonium transition into gg-pair are accounted for by normalising the short-distance 
part of the hybrid light-cone wavefunctions to the width of heavy quarkonium decay 
into leptons. The analysis of [5] indicated that once the short-distance corrections were 
included (for transverse distances less than b — 0.1 GeV -1 ) several potential models 
produced similar results. As a result, we expect that the resultant model uncertainty 
is less than that from other sources and only consider the potential of [8]. 
In the photoproduction limit we have 

<+>" - vp) - ISwS c « ! = °» • < 2 - 3 > 

where, as usual, an exponential decay in \t\ is assumed, with B D V the slope parameter. 
Since these parameters have not yet been measured, we are forced to estimate them. 
The established trend is that the slope parameter decreases with increasing vector me- 
son mass. This is a geometrical effect: most of Bp comes from the proton end of the 
ladder where the typical transverse momenta are much smaller, the upper end con- 
tributes progressively less as the typical transverse momentum scales involved increase. 
On this basis we expect B Dr to be a little less than that measured in photoproduction 
of J/ip: B Dy j/^, = AA ± 0.3 GeV~ 2 (HI). For the purpose of this paper we assume, 
for each state, the same value as that measured in electroproduction of J/ip by HI, 
Bd,v — 3-9 GeV -2 (see [9] and references therein), since the typical transverse scales 
at the upper vertices in these processes are similar. Of course, this introduces a further 
overall uncertainty in the normalisation of the cross sections of about 10 — 20%. 

3. Rescaling: calculation of Q 2 e fj 

The procedure that we use for determining the scale, Qlff, in eq.(2.3) is similar to that 
described in sec. (Ilia) of [4]. The amplitude the for photoproduction of T in co-ordinate 



space is given by: 
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where, 

a(b 2 ) = ^-b 2 a s (b 2 )xg(x,b 2 ) (3.2) 

is the cross section for scattering a qq dipole of transverse size b 2 off the proton (bt is 
the conjugate variable to k^). This quantity is expected to be universal, ^^(z, bf) — 
Ko(bm run ) and iPv,t{z, h) are the light-cone wavefunctions of the transversely-polarised 
real photon and heavy vector meson in configuration space, respectively. 

In order to proceed we need to establish the relation between transverse sizes and 
momentum scales, such as Qlff, to determine the scale at which we must sample the 
gluon density and a s in eq.(3.2). We establish this relation using the expression for the 
longitudinal structure function written in 6 t -space: 

/poo 
dz J d 2 b t a(b 2 )\^, L (z,b t )\ 2 . (3.3) 

We now assume that for a particular fixed x, Q 2 the product of transverse size and 
the momentum scale squared is constant, bf Q 2 = A, with b\ measured in GeV -2 = bf 
(fm 2 )/(hc) 2 , with he = 0.197 GeV fm. We establish this constant from the typical or 
average transverse sizes contributing to the Vintegral above. The average b\ is defined 
to be that value up to which we must integrate in order to reach half of the full b 2 t 
integral in eq.(3.3). This average value is fed back into the integral via A =<b 2 > Q 2 . 
The integral is recalculated, sampling the gluon and a s at A/ <b 2 > and a new median 
bf is found. The procedure is then iterated to convergence. It turns out that the 
resultant A depends weakly on x,Q 2 . Table. (1) shows A for a wide range of these 
kinematic variables (in [5] a constant value of A = 8.5 at x = 1CT 3 is used for J/ip). 

Now we have established the relation between transverse size and momentum scales 
we can use it to set the scale of <j{b 2 t ) in the 6 r integral in eq.(3.1). This allows us to 
establish the effective scale for the production of each state, by a procedure similar to 
that used above. With some starting assumption for Q 2 jf we use A(x, Q 2 fj)/b 2 as the 
scale at which xg, a s , m 2 un are sampled in the integral. This requires some regulation at 
very small and very large values of bf, however, the contributions to the overall integral 
from these regions are negligible. We establish the median b 2 =<b 2 > for this integral, 
as above. This is then fed back in via Q 2 fj{new) = \{x,Q 2 ^{old))/ < b 2 > and the 
procedure is iterated to convergence. 

In this way the value of the effective scale depends on the dominant values of b\. 
Dipoles in the photon which most closely correspond to the relevant s-wave state, (have 
the "right transverse size" and momentum sharing z) contribute most to the production. 
The rescaling procedure is designed to reflect this. The scale varies depending on 
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Table 1: \(x,Q 2 ) as a function of x and Q 2 (GeV 2 ). The upper values refer to CTEQ4L 
partons and the lower to MRSTLO. 
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Table 2: The effective scale, Qlff in GeV 2 , for each s-wave T-state as a function of x 

the state in question since the relevant light-cone wavefunctions weight the integral 
in eq.(3.1) differently. The resulting effective scale, Q\$p depends weakly on x and 
strongly on the state concerned (see table. (2)). 

The fact that Q^ff increases with the mass of the state agrees with [4, 5] and is 
perhaps counter-intuitive given that the quarks are less tightly bound in the higher 
s-wave states, leading to larger typical sizes. It is worth remembering, however, that 
the precise momentum scale is governed by a convolution of the photon and vector 
meson light-cone wavefunctions, and that higher states also contain nodes which further 
confuses the issue from an intuitive point of view. In fact, the typical transverse 
momenta contributing to the production of these higher mass states is larger (see 
table. (3)) and this is reflected in the larger Qlff- The precise values of this effective scale 
are sensitive to the details of the vector meson wavefunction. The outlined rescaling 
procedure represents a reasonable estimate given the current knowledge of this quantity. 



4. Calculation of the effect of skewedness 



Recently, there has been considerable progress in calculating and understanding skewed 
parton distributions, which probe new non-perturbative information about hadrons 
and are a generalisation of conventional parton distributions (for an extensive list of 
references see [10]). The former replace the latter in expressions for hard exclusive 
cross-sections. They are relevant for a wide range of hard exclusive processes such 
as deeply-virtual Compton scattering, diffractive dijets in photoproduction, and most 
importantly for our purposes, photo- and electroproduction of vector mesons. 

In the latter processes, the skewedness, S, arises from the need to convert a space-like 
q 2 = —Q 2 , into a time-like M v , and is given by the difference in momentum fractions 
carried by the outgoing (xi) and returning (x 2 ) gluons (see fig.(l)), 

Xl W 2 + Q 2,X2 W 2 + Q 2 ' 1 } 



M 2 +Q 2 

where M 2 q is the mass of the intermediate gg-state. In terms of the light-cone variables 
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For photoproduction of T at HERA, 5 = M\jW 2 is small and lies in the range 
{0.0011, 0.017}, although it is an order of magnitude larger than for J/^-photoproduction 
in the same energy range. 

It was argued in [11] that for small x and for Qq, where parton densities weakly 
depend on x, the dependence on 5 can be neglected. It was further demonstrated 
in [11] that for large Q 2 and small x the main contribution to the skewed parton 
densities comes from the parton densities at Ql scale where x ^> x, for which the 
skewedness is very small. Hence the answer for these kinematics does not depend on 
this approximation. Freund and Guzey [12] have modified the DGLAP-evolution [13] 
package of the CTEQ collaboration [14] , which is based on a numerical grid integration, 
to produce the skewed gluon functions, G$(xi,Q 2 ), for any X\ at a fixed value of 5 (it 
is straightforward to write an interpolating subroutine to obtain G for any 5). We use 
this code and the approximation that the skewed and conventional distributions are 
the same at the starting scale, Q^. 

It is then sufficient to replace the conventional splitting functions with their skewed 
generalisations in the evolution (i.e. P a b(x) — > P a b(^,x)). Unfortunately the latter 
are only known to leading-log accuracy at present (although, very recently progress 
has been reported [15] on the next-to-leading order). For consistency one is forced 



to use only leading-log parton distributions at the input scale; we use the two most 
recent leading-order distributions, i.e. CTEQ4L [16] and MRSTLO 1 [17], and evolve to 
leading-log(Q 2 ) accuracy. The code produces the skewed gluon distribution G$(xi, Qeff) 
to replace the ordinary gluon density, xg(x, Q\ff) of eq.(2.3). The two distributions are 
shown in fig. (2), in the relevant x-range for Q 2 e ff = 40 GeV 2 characteristic of T(ls)- 
production (see table. (2)). This replacement leads to an overall enhancement factor of 
about (1.6) 2 ~ 2.6 for T(ls,2s,3s) cross-sections. 
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Figure 2: A comparison of skewed and conventional gluon densities at a scale of Q 2 = 40 
GeV 2 , the effective scale for Y-photoproduction. Conventional leading-order input distribu- 
tions are used in both cases (at starting scales Qq = 1.6 and 1.0 GeV 2 for CTEQ4L [16] and 
MRSTLO [17], respectively). The skewed distributions are the upper curves and have fixed 
Xl /S = 1.2. 

It is necessary to determine the appropriate value of x\ to use. We estimate this by 
calculating the average values of transverse momentum and light-cone momentum shar- 
ing <k t > and < z > contributing to M 2 ^ in the relevant loop integral. A more accurate 
calculation would integrate explicitly over M 2 ^, but in this case all direct connection to 
the measured forward distribution, and hence predictability, would be lost (until the 
input densities for skewed evolution have been explicitly measured themselves). 

The amplitude of eq.(3.1) is used to calculate the relevant value of M 2 ^ for the 
particular value of W. First of all we calculate <z>. We then integrate over z, take 
a Fourier transform back to fc t -space and calculate the median value of the /^-integral 
by the method outlined above for calculating medians of integrals. Both < z > and 
< A; 2 > are then fed into eq.(4.3) to determine M 2 ^. The range of values of < kf > and 
typical <z> for each s-wave state are shown in table. (3). The result depends weakly 
on energy, W. Table. (4) shows the resultant values of Xi/8 — M^/M v as a function of 

1 We thank R. Roberts for providing the parameters at the starting scale. 



W for the three s-wave states. The code for calculating the cross sections automatically 
takes this variation into account. 

Note that as soon as x±/5 is small we are not sensitive to a specific value of this 
ratio since in this limit the ratio of skewed and diagonal densities weakly depends on 
Xi/5. This justifies the use of the conventional distributions in a when calculating M 2 ^. 
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Table 3: Averages of k 2 and z for each s-wave T-state. 
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Table 4: The ratio x\/8 for each s-wave T-state as a function of W . 



5. Calculating the real part of the amplitude 

It is convenient to use the dispersion representation in energy to calculate the real part 
of the scattering amplitude. In doing so we automatically take into account a possible 
contribution of the region of X\ > 0, x 2 < in the real part. Previously, the ratio of real 
to imaginary parts of the amplitude, f3, was calculated using the approximate solution 
of this dispersion relation: 

_ ir dln(xg(x,Q 2 )) 
* "2 dln(l/x) ' [ ' 

first derived in [18], and used in [4-6]. This is appropriate at small x and fairly low 
scales where the gluon density and hence the imaginary part of the amplitude may be 
approximated by a single power in energy (ImA oc xg/x oc (W 2 ) a ). The result is then 
P = | (a — 1), for a sufficiently close to 1. 
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Table 5: Relative contribution of the real part of the amplitude, /3 2 , as a function of W for 
photoproduction of T(ls), T(2s), T(3s). For comparison the last column shows the values of 
(3 2 obtained for T using eq.(5.1), sampling the CTEQ4L gluon at = 40GeV 2 . 

At the larger values of x and higher scales, Q\fp relevant to photoproduction of T- 
states it is necessary to include an additional sub- leading power in 1/x or equivalently 
in W 2 . For fixed X\/8 (we use the values at W = 180 GeV, see table. (4)) and Qlff the 
skewed gluon distribution is a function of W 2 only. We perform a two-power Regge- 
type fit to this distribution, at the relevant fixed effective scale, over the whole range 
in W 2 using MINUIT [19]: 

G s (x 1 , Q 2 eff ) = aW 2b + cW 2d . (5.2) 

We may then use analyticity of the amplitude directly to construct a dispersion 
relation and hence determine ft, 



ReA aW 2 ^ cot(^) + cW 2 ^ cot(^) 
f[ ] ~ ImA~ aW 2(b + i) + cW 2{d+i) ■ ^-3) 

Table. (5) shows how (3 2 chang function of W, these values are used in the cross 

sections of eq.(2.3). For a comparison, we also show the values of f3 2 for T using eq.(5.1). 
As expected these differ most for smaller values of W . 



6. Results 



The cross section for photoproduction of the T states is 

3n 3 T v M^(l + p 2 ) 
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The effective scale, Qljf, x±/S and (3 2 all depend on the state concerned and weakly 
on energy as indicated in the tables. (2,4,5) above. The overall suppression factor, 
C(Q 2 = 0) is given by 
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The resolution of the muon chambers at the HERA experiments is such that the 
three contributing s-wave states T,T',T" are not resolved [7]. Hence, what is actu- 
ally measured is the sum of the production rates of the mesons times their respective 
branching fractions to muons (Br(T (ns) — > = 2.48, 1.31, 1.81 % for n = 1, 2, 3, 

respectively [21]). In fig. (3) we present our predictions for this measured quantity using 
MRST (upper solid curve) and CTEQ (lower solid curve) input parton distributions. 
This illustrates the agreement of the model with the data collected so far. Also shown 
(dashed curves) are the predictions for T(ls) alone for the two input distributions. 
There are two reasons why the MRST curves should be higher than CTEQ curves. 
Firstly, the skewed parton distributions are higher (see Fig. (2)). Secondly, smaller val- 
ues for Aqcd are used in MRSTLO than in CTEQ4L, this makes the value of a s larger 
at a given effective scale. 
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Figure 3: Sum of cross sections times branching ratios to muons for T(ls), T (2s), T (3s) 
as a function of energy. The dashed curves are for T(ls) alone. The upper (lower) curves in 
each case correspond to MRSTLO (CTEQ4L) partons at the starting scale. Also shown are 
the data from ZEUS and HI (preliminary), at their respective mean energies, with systematic 
and statistical errors added in quadrature. 



Fig. (4) shows our prediction for the photoproduction of T(ls) as function of energy 
in the HERA range, using both input distributions. A very steep rise is expected, as 
can been seen from the figure, this corresponds approximately to 1U 1 ' 7 over the range 
shown, i.e. almost a full power in W stronger than that seen in J/ip production (a 
recent fit [9] to ZEUS and HI data revealed a power of 0.8 ± 0.1). This very steep rise 
is due to the sampling of the gluon at the large scale, Q 2 ef j ~ 40 GeV 2 , where it is 
rising steeply with energy. In practice the steepness of the observed rise may prove to 
be a useful way to discriminate between models which have rescaling and those that 
do not, since for a fixed range in x the steepness increases with scale. 
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Figure 4: Cross section for photoproduction of T(ls) in the HERA energy range, compared 
to values quoted by the HERA experiments [7] at their respective mean energies. An addi- 
tional curve, AW 1 - 7 , with A normalised to the CTEQ4L cross section at W = 80 GeV, is also 
shown to indicate the very steep rise with energy. 
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Figure 5: Ratios of cross-sections times branching ratio to muons: Rz s 
R 3s = aB(3s)/aB(ls). 



aB(2s)/aB(ls), 



Figure. (5) shows predictions for ratios of the products of cross-sections and branch- 
ing ratios to muons, for the T', T" relative to T as a function of energy. The different k t - 
suppression factors for the different states (for CTEQ4L C{Q 2 = 0) = 0.290, 0.139, 0.0820, 
for MRSTLO C(Q 2 = 0) = 0.283,0.138,0.0817, respectively) imply that the relative 
rates of the three states are not expected to be the same as those found at CDF [20], 
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and used by ZEUS and HI to unfold and produce a cross section for T(ls) production. 
The CDF data implies that T(ls) is responsible for about 70% of the signal whereas 
our calculation indicates a larger value of about 85% (see also fig. (3)). For this reason 
a comparison of the predictions in Fig. (4) with the values quoted by the experiments 
(also shown), must be done with care. 

Although the states cannot be separated this predominance of T(ls) should reflect 
itself in the average observed value of the mean mass of the muon pair in the signal. 

7. Conclusions and discussion 

We have presented predictions for the cross sections of diffractive photoproduction of 
T-states in the HERA range, based on a previous analysis [5]. In this paper two new 
effects, the off-diagonal nature of the amplitude and the largeness of it's real part are 
found to be important and lead to a significant enhancement of the normalisation of 
the cross sections. We perform a reanalysis of rescaling procedure used to set the scale 
of the production of the three s-wave states. We find good agreement with the first 
data from HERA. One striking prediction of the analysis is a very strong rise of the 
cross sections with energy which is driven by the skewed gluon density (see eq.(6.1) 
and fig. (2)) 

Throughout the paper we have used the hybrid light-cone wavefunctions suggested 
in [5], which were constructed to agree with quarkonium wavefunctions at large dis- 
tances and QCD at small distances. These wavefunctions have significant Fermi motion 
of the quarks inside the bound states. Gauge invariance demands that if one has finite- 
k t one also needs gluonic degrees of freedom. So far only the contribution to the cross 
section of the lowest order Fock states of the vector meson and photon (\bb>) have been 
considered. The next step is to calculate higher order Fock states which also include 
the gluon degrees of freedom. The connection between these light-cone wavefunctions 
and quarkonium wavefunctions derived from solving the Schrodinger equation for a 
particular potential remains an open question, which such a calculation would begin to 
address. 
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